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In order to illustrate the capabilities of COMPASS using a hadronic beam, I review some of the 
azimuthal asymmetries in hadronic collisions, that allow for the extraction of transversity, Sivers 
(3JT)j and Boer-Mulders functions, necessary to explore the partonic spin structure of the nucleon. I also 

^ . report on some Monte Carlo simulations of such asymmetries for the production of Drell-Yan lepton 

' pairs from the collision of high-energy pions on a transversely polarized proton target. 



*5 ' !• INTRODUCTION 

Oh 1 The recent measurement of Single-Spin Asymmetries (SSA) in semi-inclusive lp' — ► I'irX Deep-Inelastic Scattering 
(SIDIS) on transversely polarized hadronic targets [l|, SHI, El, has renewed the interest in the problem of describing 
i— —i[ the spin structure of hadrons within Quantum Chromo-Dynamics (QCD) Experimental evidence of large SSA in 
hadron-hadron collisions was well known since many years [f| @] , but it has never been consistently explained in the 
context of perturbative QCD in the collinear massless approximation Q . Similarly, in a series of experiments on high- 
energy collisions of pions and antiprotons with various unpolarized nuclei [f| [T3, EH E3] the cross section for Drell-Yan 
events shows an unexpected largely asymmetric azimuthal distribution of the final lepton pair with respect to the 
, production plane, which leads to a bad violation of the socalled Lam- Tung sum rule Again, in the context of 
collinear perturbative QCD there are no calculations able to consistently justify this set of measurements [lil fTBl [l6| ■ 
The idea of going beyond the collinear approximation opened new perspectives about the possibility of explaining 
these SSA in terms of intrinsic transverse motion of partons inside hadrons, and of correlations between such intrinsic 
\ transverse momenta and transverse spin degrees of freedom. The most popular examples are the Sivers [l7j and the 
Collins [l8| effects. In the former case, an asymmetric azimuthal distribution of detected hadrons (with respect to 
the normal to the production plane) is obtained from the nonpcrturbative correlation px x P ■ St, where px is the 
intrinsic transverse momentum of an unpolarized parton inside a target hadron with momentum P and transverse 
polarization St- In the latter case, the asymmetry is obtained from the correlation k x PhT • St, where a parton 



> 
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x. 

with momentum k and transverse polarization St fragments into an unpolarized hadron with transverse momentum 
PhT- In both cases, the sizes of the effects are represented by new unintegrated, or Transverse-Momentum Dependent 
(TMD), partonic functions, the socalled Sivers and Collins functions, respectively. 

However, SSA data in hadronic collisions have been collected so far typically for semi-inclusive pp^> — > h^X 
processes, where the factorization proof is still under debate (see Ref. [19( and references therein). On the contrary, 
the theoretical situation of the SIDIS measurements is more transparent. On the basis of a suitable factorization 
theorem [2(| Hl[ , the cross section at leading twist contains convolutions involving separately the Sivers and Collins 
functions with different azimuthal dependences, sin(</> — 4>s) and sin(</> + <ps), respectively, where 4>,4>S, are the 
azimuthal angles of the produced hadron and of the target polarization with respect to the axis defined by the virtual 
photon [22| ■ According to the extracted azimuthal dependence, the measured SSA can then be clearly related to one 
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effect or the other [l|, H| . 

Similarly, in the Drell-Yan process HiH\ — > l + l~ X the cross section displays at leading twist two terms weighted 
by sin(0 — <f>g) and sin(0 + 0s), where now 0, 05, are the azimuthal orientations of the final lepton plane and of 
the hadron polarization with respect to the reaction plane [23j]. Adopting the notations recommended in Ref. [HI, 
the first one involves the convolution of the Sivers function f^p with the standard unpolarized parton distribution 
fi. The second one involves the yet poorly known transversity distribution hi and the Boer-Mulders function hj~, 
a TMD distribution which is most likely responsible for the violation of the previously mentioned Lam- Tung sum 
rule (23j . Hence, a simultaneous measurement of unpolarized and single-polarized Drell-Yan cross sections would 
allow to extract all the above partonic densities from data [25|, Hfjj]. Both hi and describe the distribution of 
transversely polarized partons; but the former applies to transversely polarized parent hadrons, while the latter to 
unpolarized ones. On an equal footing, f^j, and fi describe distributions of unpolarized partons. The correlation 
between pt and St inside f^~ T is intuitively possible only for a nonvanishing orbital angular momentum of partons. 
Then, extraction of Sivers function from SIDIS and Drell-Yan data would allow to study the orbital motion and 
the spatial distribution of hidden confined partons (27[, as well as to test its predicted sign change from SIDIS to 
Drell-Yan, a peculiar universality feature deduced on very general arguments that holds also for hi |28|. On the 
same footing, the prediction that the first moment of transversity, the nucleon tensor charge, is much larger than its 
helicity, as it emerges from lattice studies [H,|33|, represents a basic test of QCD in the nonperturbative domain. 

In a series of papers [25|, l3ll . l32l . |33| , numerical simulations of SSA in Drell-Yan processes were performed using 
transversely polarized proton targets, and (anti)proton and pion beams, in different kinematic conditions correspond- 
ing to the setup of GSI, RHIC, and COMPASS. After briefly reviewing the formalism in Sec. [Til and the strategy 
for Monte Carlo simulation in Sec. IIIH in Sec. IIVI we will reconsider the sin(0 ± </>g) SSA for the n^p' _> ^+ p- x 
process using the pion beam at COMPASS, in order to verify if the foreseen statistical accuracy is enough to test the 
predicted sign change of and to extract information on /if. 

In the SIDIS production of pions on transversely polarized targets, the above mentioned Collins effect represents the 
most popular technique to extract the transversity hi from a SSA measurement. However, it requires the cross section 
to depend explicitly upon PhT, the transverse momentum of the detected pion with respect to the photon axis 22]. 
This fact brings in several complications, including the possible overlap of the Collins effect with other competing 
mechanisms and more complicated factorization proofs and evolution equations 0, [2l[ . Semi- inclusive production of 
two hadrons [34], [35[ offers an alternative framework, where the chiral-odd partner of transversity is represented by the 
Dihadron Fragmentation Function (DiFF) (36[ , which relates the transverse spin of the quark to the azimuthal 
orientation of the two-hadron plane. This function is at present unknown. Very recently, the HERMES collaboration 
has reported measurements of the asymmetry containing the product hi [37j • The COMPASS collaboration 
has also presented analogous preliminary results |38[. In the meanwhile, the BELLE collaboration is planning to 
measure in the near future [3^|. A spectator model calculation of leading-twist DiFF has been built by tunin, 



the parameters on the output of pion pair distributions from PYTHIA, adapted to the HERMES kinematics 4C]. 
Evolution equations for DiFF have been also studied, including the full dependence upon the pair invariant mass 41 1. 

In Sec. [V] we will reconsider the proposal formulated for the first time in Ref. [42j , namely to use hadronic collisions 
on a transversely polarized proton target and inclusively detect one (or two) pairs of pions. By combining unpolarized 
and polarized cross sections, in principle it is possible to "self-determine" all unknown DiFF in one experiment and 
to extract the transversity hi in a clean way. We think that the possibility should be explored to perform such a 
program at COMPASS using a pion beam. 




FIG. 1: The Collins-Soper frame. 
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II. THE FORMALISM FOR THE DRELL-YAN CROSS SECTION 



In the following, we will just sketch the kinematics and the necessary formulae for describing the H\H\ — ► l + l~ X 
cross section at leading twist. We forward the interested reader to Ref. [33] for the technical details. 

In the collision of two hadrons with momentum P 1; P 2 , and masses Mx,M 2 , the center-of-mass (cm.) square energy 
s = (Pi + P 2 ) 2 is available to produce two leptons with momenta k\, k 2 , and invariant mass q 2 = M 2 = (ki + k 2 ) 2 ■ In 
the limit M 2 , s — » 00, while keeping the ratio < t = M 2 / s < 1 limited, a factorization theorem can be proven [43[ 
ensuring that the elementary mechanism proceeds through the annihilation of a parton and an antiparton into a virtual 
photon with time-like momentum q 2 . The partons are characterized by fractional longitudinal momenta xi,x 2 , and 
spatial trasverse components PiT with respect to the direction defined by Pi (i — 1,2), respectively. Momentum 
conservation implies qt = Pit +P2T- If Qt 7^ the annihilation direction is not known. Hence, it is convenient to 
select the socalled Collins-Soper frame [4J] described in Fig. [TJ The final lepton pair is detected in the solid angle 
(9, </>), where, in particular, (j) (and all other azimuthal angles) is measured in a plane perpendicular to the indicated 
lepton plane but containing h = qr/lqrl- 

The full expression of the leading-twist differential cross section can be written as [23| 
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where a is the fine structure constant, dQ, — sin 
angle of the transverse spin of hadron i, and 



, e q is the charge of the parton with flavor q, ^>Si is the azimuthal 
\-y + V 2 ) C - J(l + cos 2 ^) B(y)=y(l-y) ^ \ sin 2 9 . (2) 



The convolutions arc defined as 



T 



DF^DFUhP] 



(X) 

dp 1T dp 2T 5 (pit + Pit ~ qr) DF 1 (x 1 ,p 1T ;q/H 1 ) DF 2 (x 2 ,p 2T ;q/H^') + (q <-> q) 



Parton distribution functions can be obtained from the color-gauge invariant quark-quark correlator 



Hx, s) = / ^- e ™ p+ z- (p, s\m u [0A -] m\p, S) 



(3) 



(4) 



J+=£ T =o 



where 



U %t] =Ve- i9 ti dw - A ^ (5) 

is the socalled gauge link operator, connecting the two different space-time points 0,£, by all the possible ordered 
paths ("P[..]) followed by the gluon field A, which couples to the quark field ip through the constant g. Only A + 



4 



p" 












V 1 










p" 



FIG. 2: A possible representation of the color-gauge invariant quark-quark correlator. 



gluons appear at any power in the expansion of the exponential, leading to an intuitive representation of the gauge 
link operator as a sort of residual interaction between the parton and the residual hadron (see Fig. [21). When also 
including the pt dependence in Eq. UJ the integration path of the gauge link is more complicated and involves At 
gluons at light-cone infinity [45| . For the specific case of interest of hadrons with transverse polarization St, the 
quark-quark correlator can be parametrized as 



1 f P v p p S C! p p P v 

<f>(x,p T ,S T ) = {/lfoPr) If + AtO^Pt) WaT*" Mg T + K (x,p T ) o> 

-hxix^^ia^S^P" ^ - hi T (x,p T ) ~ JJT S T P V ^ io^l^ , (6) 



where Mh is the parent hadron mass. 

The TMD functions of Eq. ([TJ) can thus be obtained by suitable projections 



^ (x, p T ,S T ) = \J dp- Tr [r §(p, P, S T )} 



(7) 



p+=xP+ 



In fact, following Ref. [241 ] we can define the number density f qST /HS T °f partons with flavor q and transverse 
polarization st in a hadron H with transverse polarization St- 



JqST/HST r, l T ~ w 



1 I ,«/ \ Pxpt-St ,±q, Pxpt-st 
4 <j fi(x,PT) ~ fi T ( x >PT) ^- v (2:,Pt) — 

,2 



, uii \a a?/ (pt-StPt-st P t a \\ (Q , 
+h q 1 (x,pT)S T -ST + h 1T 1 (x,p T ) y-^j^ — m| t ' St Jj ' ® 

where n_ = (1,0, Ot)- It is easy to verify that the unpolarized parton distribution is obtained by summing the 
number densities for all possible four combinations of transverse polarizations: ff = / ? | /h]+ fqi/H]+ fq\ /HJ. + fql/Hi- 
Similarly, the Sivers function emerges as the unbalance between the probabilities of finding unpolarized partons inside 
hadrons transversely polarized in opposite directions, namely 

(/Vr/pT + fql/pl) ~ (/aTM + fiihi) = Ulv\ ~ UlvX = ~ f it ( x >Pt) jj^ — - • (9) 

The Boer-Mulders function and the transversity emerge from the other two combinations 

_i_ P x pt ■ st 

(/•jT/pT + fql/pl) ~ (fql/pl + fql/pl) = Uxiv ~ fql/p = -h^{x,p T ) jj- , (10) 

(/ 9 TM _ ftf/Pl) - (fqlM ~ Ulfr\) = h l( x > Pt) S t • s T + h^(x, Pt) -JfT- ST ~ JfT S t ' s T^j -(H) 
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III. MONTE CARLO GENERATION OF DRELL-YAN EVENTS 

The implementation of Eq. ([1]) in the Monte Carlo generator is represented as [25[ : 

da 



dVtdx\dx 2 dqT 



= K-\T{q T , Xl ,X2,M)\ 2 c t (q T ,x 1 ,x 2 )S l {e^,cj )s2 ) , (12) 



where the event distribution is driven by the elementary unpolarized annihilation |T| 2 . We assume a factorized 
transverse- momentum dependence in each parton distribution such as to break the convolution T in Eq.Q, leading 
to 

|T| 2 « A(q T , X U X2,M) F{x u x 2 ) , (13) 

where = \qr\- The function A is parametrized and normalized as in Ref. [111 ], where high-energy Drell-Yan -k — p 
collisions were considered. The average transverse momentum turns out to be (qr) > 1 GeV/c, which effectively 
reproduces the influence of sizable QCD corrections beyond the parton model picture of Eq. fl}. It is well known [46[ 
that such corrections induce also large K factors and an M scale dependence in parton distributions, determining their 
evolution. As in our previous works [25|, [3l|, HH, [33| , we conventionally assume in Eq. ([12")) that K = 2.5, but we stress 
that in an azimuthal asymmetry the corrections to the cross sections in the numerator and in the denominator should 
approximately compensate each other, as it turns out to happen at RHIC [13] and JPARC [48| cm. energies. Since the 
range of M values here explored is close to the one of Ref. [11(, where the parametrization of A, F, and c, in Eq. (JT3J) 
was deduced assuming M-independent parton distributions, we keep our same previous approach (25l. l3ll l32l. l33j] and 
use 

2 

F(x 1 ,x 2 ) = ^— ^ e 2 flix^q/Hx) f?(x 2 ;q/H 2 ) + (q <-> q) , (14) 

9 

where the unpolarized distribution f1(x) for various flavors q = u,d, s, is taken again from Ref. [ll| . unless an explicit 
different choice is mentioned. 

The whole solid angle (8, <fi) of the final lepton pair in the Collins- Sop er frame is randomly distributed in each 
variable. The explicit form for sorting it in the Monte-Carlo is [25j, l32|, l33j 

Y,c l (qT,x 1 ,x 2 )S l (e^,cbs 2 ) = (l + cos 2 g)+ KgllfglgO sin 2g CQs20 

i=l 

+ 1 S 2 t | c 3 (q T ,xx, x 2 ) S , 3 (6', 0, 0s 2 ) . (15) 

If quarks were massless, the virtual photon would be only transversely polarized and the angular dependence would 
be described by the coefficient Ci = 1, Si = 1 + cos 2 9. In the following, we discuss violations of such azimuthally 
symmetric dependence, as they appear in Eq. (|15|) . 



A. The cos 20 asymmetry 

Azimuthal cos 20 asymmetries in Eq. (fl2|) . induced by c 2 = are due to the longitudinal polarization of the 
virtual photon and to the fact that quarks have an intrinsic transverse momentum distribution, leading to the explicit 
violation of the socalled Lam- Tung sum rule [ll| . QCD corrections influence v, which in principle depends also on 



M 2 [llj . The coefficient c 2 was simulated at the GSI kinematics [25| , using the simple parametrization of Ref. [2^ 
and testing it against the previous measurement of Ref. [Io| . The pt dependence of the parton distributions in the 
nucleon was parametrized as [23| 

f!(x,p T ) = ^e-T&mx) 

TT 

h±*(x,pT) - Jf^ /ffopr), (16) 
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where ax = 1 GeV 2 , M c = 2.3 GeV, c q = 1, and Mm is the nucleon mass. Following the steps in Sec. VI of Ref. [2c 
one easily gets 



E 9 



, , h^ q (x 1 ,p 1T )h^ q (x 2 ,P2T 
2hp 1T h- p 2T - Pit • P2T ) 



+ (1 «->2) 



(# + 4M 2 ) 2 



E, ^ [/f(a:i,Pir) A 9 (asa,Par)] + (1 - 2) 

(17) 



Since the dependence of h± was fitted to the measured cos 2(f) asymmetry of the corresponding unpolarized Drcll- 
Yan cross section, which is small for 1 < qx < 3 GeV/c (see, for example, Fig. 4 in Ref. [23]), correspondingly, the 
cos 2(f> asymmetry turned out to be small for the considered statistically relevant qx range in kinematics conditions 
reachable at GSI [Hj]. 

In the literature, there are several models available for in the nucleon; in particular, in the context of the 
spectator diquark model both T-even [491 ] and T-odd parton densities I50l have been analyzed, and further work is in 
progress (5l| . Here, we reconsider the spectator diquark model of Ref. [52t| and we present some preliminary results, 
obtained in the context of an improved and enlarged approach [53| . The main features of the latter are the following. 
The gauge link appearing in the quark-quark correlator is expanded and truncated at the one-gluon exchange level 
(see Fig. [5]). The emerging interference diagrams produce the naive T-odd structures necessary to have nonvanishing 
TMD distributions like the Sivers fyj. and the Boer-Mulders hi functions. Because of the symmetry properties of 
the nucleon wave function, scalar and axial-vector diquarks are considered with different couplings to the valence 
quark left. A dipole-like form factor is attached at each nucleon-quark-diquark vertex (but other choices will be 
explored [53j]). The obtained expressions, that replace the ones in Eq. (|16p. are 

_ N 2 qS [p T + {m + M N xf] (l-xf N 2 a [p 2 T (l + x 2 ) + ( m + M N x) 2 (l-x) 2 } (1 - x) 

, = Nh M N (m + M N x) (1 - x)3 Nl M N (m + M N x) (1 - x) 2 

1 64 - 4 A%{Al+p T f 64^ Al(Al+p T f 

where 

16tt 
~3 



A 2 s/a = (1 - x)A 2 +xM 2 s/a - x(l -x)Ml , ee D = -4nC F a s » — 0.3 , (19) 



and m — 0.36 GeV is the valence (constituent) quark mass, A = 0.5 is a cutoff for high quark virtualities, Ms = 0.6 
GeV and M a = 0.8 GeV are the scalar and axial-vector diquark masses, respectively. With these choices, the 
normalizations turn out to be N^ s = 6.52 x 3/2, N% a = 30.63 x 1/2, N% s = 0, Nj a = 30.63. The sign of the obtained 
h\ q (and also of /{^?, see below) is consistent with the lattice predictions [3(| for both q = u,d. 

As a preliminary step in the analysis of azimuthal asymmetries at COMPASS, for the evolution of the parton 
densities in this context we made an educated guess, since this information is unknown for both and . We have 
integrated the px dependence away in Eq. (JTHJ) , and applied DGLAP evolution at NLO using the code of Ref. 
for unpolarized quarks and of Ref. [551 ] for transversely polarized ones. The resulting expressions have been inserted 
in Eq. (|16[) in place of the corresponding f1{x), in order to have an evolved model dependence in the behaviour in 
x, but retaining the phenomenological behaviour in px fitted to the available data. Consistently with the discussed 
approximations, also the evolution scales have been determined in an approximate way. The final scale obviously 
reflects the foreseen COMPASS setup, where Drell-Yan dileptons can be detected with invariant masses above 4 GeV; 
hence, DGLAP evolution has been stopped at Q 2 = 16 GeV 2 . The initial soft scale has been deduced by solving 
the renormalization group equations for the second Mellin moment of the non-singlet (valence) distribution, with the 
anomalous dimension consistently determined at NLO; the result is Qq 0.1 GeV 2 . In Fig. [3] the model xhj~ u (x) is 
shown at Qq (blue line) and at Q 2 (black line); the red line is a fit to the latter one with the form Nx a (l — x)' 3 , which 
is actually used in the Monte Carlo calculation of the coefficient c£ to produce the azimuthal asymmetry. 



B. The Sivers effect 



If we consider the Sivers effect in Eq. ([T|), the last term in Eq. (fT5")l becomes 

Sf v (e, (f, (f> S2 ) = (1 + cos 2 d) sin(0 - (f>s 2 ) 



(20) 



Boer-Mulders up 



FIG. 3: The spectator diquark model calculation of the Boer-Mulders function xh^ u (x) at two different scales: initial Qq —0.1 
GeV 2 (blue line), final Q 2 = 16 GeV 2 (black line). DGLAP evolution is considered at NLO level. The red line is a 3-parameter 
fit to the black line (see text). 



and the corresponding coefficient C3 reads 



cf v (qT,x 1 ,x 2 ) 



t f\{x\ , Pit) /jV^, P2T) 
HP2T ~AU 



(21) 



Eg e^[ff(xt,piT) ff( x 2,P2r) 
where the complete dependence of the involved TMD parton distributions has been made explicit. For the Sivers 



function we adopt the parametrization described in Ref. 32]. It is inspired to the one of Ref. [56j . where the 
transverse momentum of the detected pion in the SIDIS process was assumed to come entirely from the pr dependence 
of the Sivers function, and was further integrated out building the fit in terms of specific moments of the function itself. 
The x dependence of that approach is here retained, but a different flavor-dependent normalization and an explici t pj 
dependence are introduced that are bound to the shape of the recent RHIC data on pp^ — > ttX at ^/s — 200 GeV [57| , 
where large persisting asymmetries are found that could be partly due to the leading-twist Sivers mechanism. The 
expression adopted is 

fi T g (x,p T ) = N q x(l-x) M2P ° P J^ /i0, Pt) 

= N q x (1 - x) J£&* 1 e -4/<^> m , (22) 

where po = 2 GeV/c, and N u = —Nd = 0.7. The sign, positive for u quarks and negative for the d ones, already takes 
into account the predicted sign change of f^. from Drell-Yan to SIDIS [28[ . 

Following the steps described in Sec. III-2 of Ref. [32j, we can directly insert Eq. (f22|) into Eq. (f2~Tj) and get 

Siv _ f 2poq T \ 2 8N U + N d 

c 3 «x 2 (l-x 2 ) 9 ■ (23) 

The qx shape and peak position are in agreement with a similar analysis of the azimuthal asymmetry of the unpolarizcd 
Drell-Yan data (the violation of the Lam- Tung sum rule [HI]), but are induced also by the observed x F —q T correlation 
in the above mentioned RHIC data for pp^ — > ttX , when it is assumed that the SSA is entirely due to the Sivers 
mechanism. This suggests that the maximum asymmetry is reached in the upper valence region such that x F ~ X2 ~ 
(«r>/5 G3- 

C. The Boer-Mulders effect 

If we consider the Boer-Mulders effect in Eq. |T]) , the last term in Eq. (fT5"l) becomes 

5 3 BM (0, cj>, 4> S2 ) = sin 2 9 sin(0 + 0s 2 ) (24) 
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and the corresponding coefficient C3 reads 

h ± q (xi, Pit) h>l(x 2 , p 2 r) 



E 9 

cf M (<7T, Xi,X 2 ) = — 



h ■ pit 



Mi 



E g e2^[/f(a;i,p 1T )/ 1 9 (x 2 ,p 2T )] 
As in Ref . [23] , one can use Eq. (|16[) and the following parametrization for transversity, 



(25) 



h\{x,pT) = —e- a TPTh\{x), (26) 



7T 



to get 



BM 2M C 9T E q e 2 q fKx^q/Hx) h\{x 2 -q/Hl) + (q «-> g) 



3 # + 4M 2 £ ? e 2 fHx^q/H,) fl(x 2 - q/H 2 ) + (q ^ q) 

2M c g T f(x 1 ;(g)/H 1 )h 1 (x 2 ;(q)/Hl) 2M C q T fei(x 2 ; (g)/ff 2 T ) 

9 2 +4M 2 f( Xi; {q)/ Hl )f 1 (x 2 ;(q)/H 2 ) q% + 4MJ frfa; (q) /H 9 ) ' 1 j 

where the second step is justified by assuming that the contribution of each flavor to each parton distribution can be 
approximated by a corresponding average function. This approximation, together with the systematic neglect of sea 
(anti)quark contributions, has been widely explored in Ref. (58|, and it turns out to be largely justified at COMPASS 
kinematics. Because of the lack of parametrizations for h± , numerical simulations of cf M were performed in Ref. (33] ] 
by making two different guesses for the ratio hi(x 2 ; (q) ) / fx (x 2 ; {q)/H 2 ) 1 namely the ascending function ^fx-2. and 
the descending one \J\ — x 2 , that both respect the Soffer bound. The goal was to determine the minimum number 
of events (compatible with the kinematical setup and cuts) required to produce azimuthal asymmetries that could be 
clearly distinguished like the corresponding originating distributions. 

Alternatively, here we will show preliminary results for the c BM asymmetry obtained in the spectator diquark 
model [53[ by using the same strategy discussed in Sec. IIII Al Namely, we start from the ^-integrated model 
expressions for f\ and h\ q from Eq. (jT5J) , complemented by the model transversity 

hQM _ N h (m + M N xf (1 - xf N 2 a X {l- X ) 

Next, we evolve these distributions from Q\ = 0.1 GcV 2 to Q 2 = 16 GeV 2 with a (polarized) DGLAP NLO kernel. 
Then, we fit each result with different 3-parameter forms NiX at (l — x)^*-, i = 1 — 3, and we replace them in the 
corresponding ^-dependent part of the parametrizations (fl"6)) and (|26|) . The calculation of cf M then develops in the 
same way up to the first line of Eq. ([2T)l , where each x-dependent parton density is now replaced by its properly-evolved 
corresponding expression in the spectator diquark model, with no need to apply the flavor average approximation. 
Note that the ff in the numerator now does not simplify with the one in the denominator, because the former and 
the latter are replaced by the pp-integrated and evolved hf q (x,pT) and /i(£,Pt) functions of Eq. (|T5|) . respectively. 



IV. RESULTS OF DRELL-YAN MONTE CARLO SIMULATIONS 



In the Monte Carlo simulation, events are generated by Eq. (|13p . and azimuthal asymmetries are produced by 
Eq. (jTTJ) for the cos 2^ asymmetry in unpolarized Drell-Yan, by Eq. (f2"3")) for the Sivers effect, and by Eq. (|2"T]) 
for the Boer-Mulders effect, respectively, using phenomcnological parametrizations or model inputs for the involved 
(TMD) parton densities, as it is described in the previous section. The general strategy is to divide the event 
sample in two groups, one for positive values "J7" of S 2 = sin 2 cos 2cf> in the unpolarized cross section (or of 
Sf lv = (1 + cos 2 9) sin(0— 05 2 ) for Sivers effect, or Sf M = sin 2 9 sin(cj) + (j)s 2 ) for the Boer-Mulders effect), and another 
one for negative values "D" , then taking the ratio (U — D)/(U + D) [25|. Data are accumulated only in the x 2 bins of 
the polarized proton, i.e. they are summed over in the x\ bins for the beam, in the transverse momentum qT of the 
muon pair and in their zenithal orientation 9. Statistical errors for the asymmetry (U — D)/(U + D) are obtained by 
making 10 independent repetitions of the simulation for each individual case, and then calculating for each x 2 bin the 
average asymmetry value and the variance. We checked that 10 repetitions are a reasonable threshold to have stable 
numbers, since the results do not change significantly when increasing the number of repetitions beyond 6. 

The transversely polarized proton target is obtained from a NH3 molecule where each H nucleus is fully transversely 
polarized and the number of "polarized" collisions is 25% of the total number of collisions, i.e. with a dilution factor 
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TABLE I: Total absorption cross sections per nucleon and related Drell-Yan event counts per month, for pion- and antiproton- 
proton collisions at s = 200 GeV 2 and luminosity C = 4 x 10 31 cm~ 2 s~ 1 , producing Drell-Yan pairs for various invariant masses 
(see text for a discussion of the kinematics and the cutoffs). 



beam 


M (GeV) 


a tot (nb/nucleon) 


events per month 


7T 


2.5-4 (no J/i>) 


0.5 


52 000 


7T 


4-9 


0.25 


26 000 


P 


1.5-2.5 


2.4 


249 600 


P 


4-9 


0.1 


10 400 



1/4 [2^| . The beam is represented by charged pions or antiprotons such that 100 < s < 200 GeV 2 , i.e. roughly the 
cm. energy available at GSI in the socalled asymmetric collider mode [25[. The reason for employing antiproton 
beams is simply due to the lack of knowledge of in the pion. Antiproton-proton collisions involve annihilations of 
valence parton densities, as do the pion-proton ones. Hence, it is anyway useful to study the shape and size of this 
asymmetry in the same kinematic conditions, even if this kind of projectile is presently not considered at COMPASS. 
The invariant mass of the lepton pair (usually, muons) is constrained in the 1.5 < M < 2.5 GeV and 4 < M < 9 GeV 
ranges, according to the value of s, in order to explore approximately the same x valence range and to avoid overlaps 
with the resonance regions of the cc and bb quarkonium systems. 

Proper cuts are applied to the qx distribution, namely 1 < qx < 3 GeV/c, in order to filter out low-qx muon 
pairs from background processes other than the Drell-Yan mechanism, and to optimize the ratio between the absolute 
sizes of the asymmetry and the statistical errors. The resulting (qx) ~ 1.8 GeV/c is in fair agreement with the one 
experimentally explored at RHIC [57| . Whenever the 8 angular dependence is represented by a sin 2 9 function, like for 
the S2 of Eq. {TSJ) and for the Sf M of Eq. ([2"g]). the angular distribution is constrained in the range 60° < 9 < 120°, 
because outside these limits the azimuthal asymmetry is too small (25| . 

The above cuts typically produce a reduction factor of 2.5 in the initial event sample [25]. For pion beams, we 
can easily think of samples as large as 250 000 events, which are then reduced to 100 000, and further to 25 000 in 
polarized collisions by the target dilution factor. It is much more difficult to get these numbers using antiproton 
beams. However, in our simulation of pp^' — > fi + \i~ X processes we will keep them in order to make consistent 
statistical comparisons with the pion beam case. In Tab[H we list several values of the total absorption cross section 
ctot P er single nucleon, as they are deduced from our Monte Carlo for pion and antiproton beams hitting a NH3 target 
at s = 200 GeV 2 and in the above specified conditions, and producing final Drell-Yan muons with different invariant 
masses. At the luminosity £ = 4 x 10 31 cm~ 2 s~ 1 and with efficiency 1, the product <Jtot£ gives the number of Drell- 
Yan events per nucleon and per second that could be ideally reached at COMPASS. For example, with pion beams 
of 100 GeV energy hitting a transversely polarized NH3 target at s = 200 GeV 2 , and producing muon pairs with 
invariant masses in the range 4-9 GeV, it is possible to collect around 25 000 Drell-Yan events within approximately 
one month of run at the luminosity 4 x 10 31 cm _2 s _1 . 

A. The cos 20 asymmetry 

We first consider the sample of 100 000 Drell-Yan events for the pp — > p + p~X reaction at y/s ~ 14 GeV, surviving 
the cuts in the muon invariant mass, 4 < M < 9 GeV, and in the angular window, 60° < 9 < 120°. Moreover, there 
is no target dilution factor; this sample can be collected in approximately 10 months of dedicated run (see Tab. UJ). 
Events are accumulated in the X2 bins of the target proton and for each bin two groups of events are stored, one 
corresponding to positive values (U) of S 2 = sin 2 #cos2</> in Eq. lfl~5"|) , and one for negative values (D). 

In Fig. 21 the asymmetry (U — D)/(U + D) is shown for each bin x^. Average asymmetries and (statistical) error 
bars are obtained by 10 independent repetitions of the simulation. Boundary values of xi beyond 0.7 are excluded 
because of very low statistics. Blue open triangles indicate the asymmetry generated by of Eq. (|17[) . using the 
parametrization of Eq. (|16|) . constrained to reproduce the experimental data of Ref. [1 01 ] and with the further cut 
2 < qx < 3 GeV/c. Black triangles show the same asymmetry when starting from the parton densities of Eq. (fT5|) 
in the spectator diquark model, as explained in Sec. IIII Al In this case, open triangles refer to simulations with the 
same cut 2 < qx < 3 GeV/c as before, while close triangles use the wider cut 1 < qx < 3 GeV/c. 

The negligible statistical error bars are due to the large sample, which is probably unrealistic for a p beam, and 
the displayed sensitivity to kinematical cuts in qx would also not be reachable in reality. But the bulk message is 
a confirmation of the findings in Refs. [23[ and [2^]: since the px dependence of is fitted (also for the spectator 
diquark model) to the measured cos 20 asymmetry of the unpolarized Drell-Yan cross section in the experimental 
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FIG. 4: The asymmetry (U — D)/(U + D) with positive (U) and negative (D) values of cos2<^) in Eq. (|15p (see text), for 
Drell-Yan events from the pp — > [i + (i~X reaction at ,/s ~ 14 GeV, 4 < M < 9 GeV, and 60° < < 120°. Blue open triangles 
for the parametrization of Eq. (|16p , constrained to reproduce the experimental data of Ref. [13] and with 2 < qr < 3 GeV/c. 
Black open and closed triangles for the results of the spectator diquark model (|18|l with 2 < qr < 3 GeV/c and 1 < qr < 3 
GeV/c, respectively. Lines are drawn to guide the eyes. 

conditions of Ref. [lOj], which is small for 1 < qx < 3 GeV/c, the resulting asymmetry in the COMPASS conditions 
turns out unavoidably small, as it was the case for the simulation at the GSI kinematics (25| . Still, the asymmetry 
itself is nonvanishing and deserves to be measured in order to attack the problem of building a parametrization for 



In order to study the Sivers effect, we consider two samples, one of 100 000 Drell-Yan events for the 7r~pT _> ^+ ^- x 
reaction again at yfs ~ 14 GeV and for muon invariant mass in the 4 < M < 9 GeV range, and another one of 
25 000 events for the n + p^ — > /i + pT X reaction in the same kinematic conditions. Both samples can be accumulated 
approximately in the same time. As before, the transverse momentum distribution is constrained in the range 
1 < qr < 3 GeV/c, the samples are collected in x-i bins and for each bin two groups of events are stored, one 
corresponding to positive values (U) of sin(</> — <fig 2 ) in Eq. (f2T)]) . and one for negative values (D). 

In Fig. [5j the asymmetry (U — D)/(U + D) is shown for each bin X2- Average asymmetries and (statistical) error 
bars are obtained, as usual, by 10 independent repetitions of the simulation. Boundary values of X2 beyond 0.7 are 
excluded because of very low statistics. The triangles indicate the results with the n~ beam obtained by Eq. (|23|) 
assuming that changes sign from SIDIS to Drell-Yan [28[. For sake of comparison, the squares illustrate the 
opposite results that one would obtain by ignoring such prediction. Finally, the open triangles and open squares refer 
to the same situation, respectively, but for the 7T + beam. 

In the valence picture of the (7r _ )7r + — p collision where the (uu) dd annihilation dominates, the SSA for the Drell- 
Yan process induced by tt + has opposite sign with respect to 7r~ because of the opposite signs for the normalization 
Nf, f = u,d, in the parametrization (|22p . Apart for very low X2 values where the parton picture leading to Eq. ([T]) 
becomes questionable, the error bars are very small and allow for a clean reconstruction of the asymmetry shape and, 
more importantly, for a conclusive test of the predicted sign change in f^ T . 



For the Boer-Mulders effect, we have collected 200 000 Drell-Yan events for the it — > X reaction at the 

lower y/s ~ 10 GeV in order to keep a significant statistics; the muon invariant mass is restricted in the 1.5 < M < 2.5 
GeV range, with 1 < q T < 3 GeV/c and 60° < 8 < 120°. In the left panel of Fig. H the asymmetry (U-D)/(U + D) 
is shown for each bin x-i between the events with positive (U) and negative (D) values of sin(</> + 0s 2 ) in Eq. (|24"|) . 
Triangles correspond to the choice h\(x2, (q) / H^) / fi(x2, {q)/H2) = \A ~ x 2 inside Eq. (f2"7| . open triangles to the 




B. The Sivers effect 



C. 



The Boer-Mulders effect 
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FIG. 5: The asymmetry (U — D)/(U + D) with positive (U) and negative (D) values of sin(0 — (f>s 2 ) m Eq. (HDJ ( se e text), for 
Drell-Yan events from the Sivers effect in the ir^p^ — » fi + fj,~ X reaction at yfs ~ 14 GeV, 4 < M < 9 GeV, and 1 < qr < 3 
GeV/c, using the parametrization of Eq. (|22C for the Sivers function (see text). Triangles indicate 100 000 events collected in 
x 2 bins for the ir~ beam and with N u > 0; squares for N u < 0. Open triangles indicate 25 000 events for the -k + beam and 
with N u > 0; open squares for N u < 0. Lines are drawn to guide the eye. 

asymmetry 
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FIG. 6: The asymmetry (U — D)/(U + D) with positive (U) and negative (D) values of sin(<75 + </>s 2 ) m Eq. 1|24[) (see text), 
for 200 000 Drell-Yan events from the Boer-Mulders effect. In the left panel, the tt~p t — > fi + fi~ X reaction is considered at 
y/a ~ 10 GeV, 1.5 < M < 2.5 GeV, and 1 < q T < 3 GeV/c; triangles for hi {x 2 , (g}/# 2 T )/7i (x 2y {q) / H 2 ) = y/l-x 2 inside 
Eq. ([27)l (see text), open triangles for h\{x 2 , (q) / H^)/ fi(x 2 , (q)/H 2 ) = yfx^; the (q) notation represents a common average 
term that replaces each contribution in the flavor sum (for further details, see text). In the right panel, the pp^ — > fi + fi~ X 
reaction is considered at y/s ~ 14 GeV, 4 < M < 9 GeV; triangles for the results from the spectator diquark model with 
1 < qr < 3 GeV/c, open triangles for 2 < qr < 3 GeV/c. In all cases, the further cut 60° < 6 < 120° is applied. Lines are 
drawn to guide the eye. 

choice hi(i2, (q) /H^)/ fi(x%, (q)/H2) = y[x~2 (see Sec. IIII Cj) . Both choices respect the Soffer bound between h\ and 
fi. The error bars are only statistical and are small because of the copious statistics. As it is evident in the figure, 
for X2 < 0.3 it is possible to distinguish the trend of the triangles (which statistically reflects the descending trend of 
the input function y/1 — Icq) from the one of the open triangles (referred to the ascending yfx^) . We made simulations 
also at the same kinematics of the previous Sivers effect, namely for yfs ~ 14 GeV and for muon pair invariant masses 
in the 4 < M < 9 GeV range. The statistical error bars remain small and the SSA is definitely nonvanishing, but a 
clear distinction between the two trends is no longer possible. 

In the right panel of Fig. [51 the same (U — D)/(U + D) asymmetry is considered for the pp^ — > [i + X reaction 
at yfs nu 14 GeV with muon invariant mass in the 4 < M < 9 GeV range and 60° < 8 < 120°. The same sample of 
200 000 events is now produced by input parton densities from the spectator diquark model. Triangles correspond to 
the events with the cut 1 < q? < 3 GeV/c, open triangles with 2 < q? < 3 GeV/c. Again, boundary values of x-i 
beyond 0.7 are excluded because of very low statistics. Within the (statistical) error bars, there is no sensitivity to 
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the cut in the lepton pair transverse momentum. 

The overall size of the spin asymmetry is small and the reached statistical accuracy indicates that the size of the 
sample is not responsible for this feature (see also Ref. [E^|). Rather, it is the outcome of a combination of several 
sources. First of all, the phenomenological pr dependence of in Eq. (|16[) . which also affects the result with the 
spectator diquark model according to the approximation described in Sec. IIII CI it is fixed to reproduce the data of 
Ref. [loj . which indicate significant asymmetries only for very large pr beyond the range of interest here. Secondly, 
the target dilution factor: it is unavoidably introduced by the features of the actual target that effectively reproduces 
a transversely polarized proton; at COMPASS the wanted "polarized" events are 1/4 of the total number of collisions 
on the NH3 molecule. Finally, the Soffer bound: the various choices for the ratio hi(x2, (<?) / H\)/ fi( x 2> {q) /H2) on 
one side, as well as the expressions (fT8]) and (|28|) for f%,hi, and hi in the spectator diquark model on the other side, 
are tightly limited by this constraint. 



V. TWO-HADRON INCLUSIVE PRODUCTION IN HADRONIC COLLISIONS 

As anticipated in Sec. HI the most popular technique to extract the transversity hi is to measure a single-spin 
asymmetry (SSA) in SIDIS production of pions on transversely polarized proton targets. The parton fragmentation 
into the detected pion can be described by a quark-quark correlator similar to the one of Eq. © , namely 



A(z,k T ) = i^Jdk+Y, J 0^e^(O\U [ ^ ] m\Ph,X)(P h ,X\mU loM \O) 



(2*) 



k-=p-/z 



)p _ ^(z,^) If h + Hi{z > k T )a llv ^\ , (29) 

which holds, in general, for a transversely polarized parton with momentum k, fragmenting in an unpolarizcd hadron 
with momentum Ph and mass M%. Applying the same definition of projection as in Eq. ([7]), and following again the 
prescriptions of Ref. [24| . we can define the number density Di l / qST of unpolarized hadrons h in a parton with flavor 
q and transverse polarization St- 



= \ U{ Z ,k T )+H^k T ) kX ^ h ST ^ , (30) 

where n + = (0,1, Ox)- It is evident that the unpolarized "decay" function is given by the sum upon all possible 
polarizations of the fragmenting quark, namely D\ — Dn/q^ + ^h/qii while the Collins function is deduced by 

D h/ql - D h/qi = ffj q (z, kr) — . (31) 

Combining the quark-quark correlators of Eq. ([6]) and (|29p . together with the elementary cross section for the virtual 
photon absorption, it is possible to build the leading-twist expression of the SIDIS cross section for a transversely 
polarized target in a factorized picture (for the general expression, see Ref. [6fj|). The transversity hi can be extracted 
by the following SSA: 

[62) 



J d^sdPhT (d^ + da 1 ) e l fi( x ) D K Z ) 



where <fis,4>h, are the azimuthal orientations with respect to the scattering plane of the target polarization vector 
and of the hadronic plane containing P/jt, respectively. From Eq. (j32|) it is evident that this strategy, known as 
Collins effect [l8[, requires the knowledge of the entire P^t distribution of the produced hadron h, which is obviously 
not possible from the experimental point of view. But also theoretically there are some complications, since the 
factorization proof and the evolution equations of TMD parton density functions are more involved (20 . [2lj . For 
the case of hadronic collisions leading to semi-inclusive final states in specific kinematic conditions, very recentl y a n 
explicit factorization- breaking mechanism has been discussed for a fcx-dependent elementary hard cross section [19| . 
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FIG. 7: In the left panel, the graphic representation of the mixed product k x Put ■ st, which governs the Collins effect. In 
the right panel, the mixed product Pi x Pi ■ st, which describes the asymmetric distribution in of a hadron pair collinear to 
the jet axis k. 



In Fig. [3 the Collins effect produced by the mixed product k x P^t ■ st (left panel), is compared with the more 
exclusive situation where two hadrons are produced [34], |35| . In this case (right panel) , it is possible to relate the 
polarization st of the fragmenting quark to the azimuthal orientation of the plane containing the two hadrons via the 
mixed product Pi x P% ■ st- If Ph = Pi + Pi is the total momentum of the pair, and 2R = P\ — Pi the relative one, 
the asymmetry survives even if the pair is taken collinear to the fragmenting quark direction fc, or alternatively if the 
dependence upon PhT is integrated away, contrary to what happens in the Collins effect. In fact, the quark-quark 
correlator for the fragmentation of a transversely polarized quark into two unpolarized hadrons reads 



8z( 



x (P h ,R; X|^(0)W (0 ,-oo)|0) 

rr<j 



k-=p-/z 



16tt 



D 1 (z,C,k T ,R T ) + i 



' 3 j 5 Gi(z,Ck T ,R 2 T ) \ 7 +, 



t%k Ti R T 



Ml 



(33) 



where z = P h [k = z\ + Z2 is the fraction of quark energy delivered to the pair, £ = 2R / P h = [z% — z-i)jz describes 



how this fraction is split inside the pair, and ej 



The new fragmentation functions appearing in Eq. 



are named Dihadron Fragmentation Functions (DiFF) [6l|. The is the analogue of the Collins function, while 
Gj~ describes the same situation but for the fragmenting quark longitudinally polarized (a sort of analogue of the jet 
handedness, see Ref. [(32j and references therein). After integrating upon kx = —PhT/z : both disappear and we are 
left with 



1 

167T 



D^zX^Rt) 



j^U<(z,CR? r )\ 



(34) 



Apart for the usual "decay" probability D\, a potentially asymmetric term survives involving the chiral-odd DiFF 
H^, which appears as the natural partner of the transversity hi in the SSA 36] 



1 



da^ - dat \R\ E,? ?! H^ q (z, C, R 2 T ) 



sin(^ fl + 4>s) da^ 



dot x M h 



E qq elf q 1 (x)Dl(z,C,R 2 T ) 



(35) 



The equally describes the fragmentation of a transversely polarized quark into two unpolarized hadrons as H^, 
but it survives the kx integration since it is related to the asymmetry depicted in Fig. which is connected to the 
Rt vector. In this respect, any observable related to the asymmetry quantitatively described by is free of the 
problems mentioned above about the Collins mechanism, in particular about the kx factorization. 

It is possible also to expand these DiFF in the relative partial waves of the two hadrons, by suitably rewriting 



the £ dependence in their cm. frame 63j. Very recently, the HERMES collaboration reported measurements of the 
SSA in Eq. (|35[) involving only the interference between s— and p— wave components, which is necessary to generate 
naive T-odd functions like [37|]. The COMPASS collaboration also presented analogous preliminary results [HI]. 
In the meanwhile, the BELLE collaboration is planning to measure H± in the near future [39(. A spectator model 
calculation of leading-twist DiFF has been built by tuning the parameters on the output of pion pair distributions 
(proportional to D{) from PYTHIA, after adapting it to the HERMES kinematics [5(|. Evolution equations for DiFF 
have been also studied, including the full dependence upon R T , or, equivalently, upon the pair invariant mass 41 1. 

Here, we reconsider the proposal formulated for the first time in Ref. [421 ]. namely to use hadronic collisions on a 
transversely polarized proton target and inclusively detect one (or two) pairs of pions. 
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We consider first the process A + B — > (h\h2)c + X, where two protons (with momenta Pa,Pb, and spin vectors 
Sa, Sb) collide, and two unpolarized hadrons hi, hi, with momenta Pi, Pi, masses Mi, Mi, and invariant mass squared 
Mq = P c = (Pi + Pi) 2 , are inclusively detected inside the same jet. The intrinsic transverse component Per with 
respect to the jet axis is integrated over and, consequently, Pc is taken parallel to the jet axis itself. The component 
Pc± perpendicular to the beam direction (defined by Pa, see Fig. [8]) serves as the hard scale of the process; it is 
assumed to be much bigger than the masses of the colliding hadrons and of M^. The following analysis is valid only 
at leading order in l/|Pc_i_|, i- e - at leading twist. 

The azimuthal angles are defined in the hadronic center of mass as follows (see also Fig. [5]) 

, (Pa x P c ) (Pa x S a ) . , (P c xS a )-1a 

coscps A = — ■— , smcps A = — , (36) 

\P A x Pel \Pa xS A \ \Pa x P c \ \PaxS a \ 

(P B x P c ) (P B x S B ) . , (Pc x gg) • P B 
cos (f>s B = -T- — • -T- — , sm</)s B = -~- — — ^ z-t> ( 37 ) 
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x S A \ 


(Pc x S B ) ■ 


Pb 


\P B x Pel -Pb 


x S B \ 


(Pa x R c ) 


Pc 


\PcxP A \\Pc 


x R c 



cos (pR — — • — , sm0i? = — - . (38 

' \P c xP A \ \P c xR c \ \P C x P A \\Pc x Rc\ 

The partons involved in the elementary scattering have momenta p a — x u Pa, Pb = x^Pb, and p c — Pc / z c . Also 
the usual Mandelstam variables have a counterpart at the elementary partonic level: 

s~x a XbS, t=—t, u = —u. (39) 

Zc Z c 

Conservation of momentum at the partonic level implies that 

-xc^t^^ x( , x a t + x b u \ CM ( \P C ±\ x a e-*a+x h e*° \ _ 

S d(S + t + U) = Z c Z c -\ = Z c Z c 7= = Z c (Zc - Z c ) , (40) 

\ X a X b S J V V s x aXb J 

where r\c is the pseudorapidity with respect to P4. 
Therefore, the cross section can be written as 



dcr „ , „ , v-^ 1 \ ^ f dx a dxb _ , „ . 

a,b,c,d (all x's) (41) 
16nS 2 M X c ,X d ;X*,X b M x' a ,xWaX' d Ac ( fc ' R-CTi 4>Rc)x' c x c 6 x' d x d . 



c%7 d|Pc_i_ I rfCc dPcr <#s A d(j> Sl 

1 



where the indices x's refer to the chiralities/helicities of the partons, and the partonic hard amplitudes M(s,i,u) can 
be taken from Refs. [42, 64]. In Eq. (|4ip. $ and A are the same parton-parton correlators of Eqs. ^ and (|33p . after 
integrating upon the intrinsic tranverse momenta and using the parton helicity basis representation. 
For initial quarks with flavor a we have [65l . |66| 

^ , a . ( fx(x a )+S AL gi(x a ) \S AT \e^Ah x {x a ) \ 

<5>a(Xa,S A ) X > aXa = I j • ( 42 ) 

V S AT e^Ahi(x a ) h(x a )-S AL gi(x a ) 
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The correlator is written in the helicity basis where Pa defines the z axis, with Sal and |SUt| indicating the parallel 
and transverse components of the polarization Sa with respect to Pa, and the x axis is oriented along Pc±- Similarly 
for quark b in the hadron B. When the initial partons are gluons, in a spin-^ target the <!> is diagonal because of 
angular momentum conservation [67t |68| . 

If the final fragmenting parton c is a quark, we have 

1 / D 1 ie^Ko %2ll H < \ 

where the helicity basis is now choosen with the z axis along Pc and the a: axis along the component of Pa orthogonal 
to Pc (see Fig. [5]). For final fragmenting gluons, because of angular momentum conservation the A c can be obtained 
in analogy to the case of gluon distributions in spin-1 targets [69(: 

i / G ie 2l ^c J5^lJ! 5G< \ 

4tt ^ _ ie -2^ Kc 5G < G ) 

The functions G and <5G^ describe the decay into two unpolarized hadrons of an unpolarized and a transversely 
polarized gluon, respectively, where by "transverse polarization" we mean linear polarization along two independent 
directions transverse to z, as in the quark case (42j . Unfortunately, 5G^ cannot appear in connection with the quark 
transversity h\ because of the mismatch in the units of helicity flip between a spin-i and a spin-1 objects, leading 
to the exp(z<^R c ) and exp(2i()> Rc ) dependences in Eq. (|4*2"|) and Eq. (14~4"|) . respectively. It can be coupled to the gluon 
transversity 5G, but only in a target with spin greater than i |69l |. 

When considering the process pp^ — ► (h 1 h 2 )X, the most interesting SSA is 

An^cAPc^XcRItArc^Ss) = daT - d(ji = , (45) 



where 



dauu = 2\P C1 _\ £ l d ^n{^)St{x b ) d ^^D 1 {z c ^Rl T ), 

da UT = 2|P C J ^ |5 BT |sin(^ B -^ c ) / ^p- f?(x a ) h\{x b ) dA& ^ d H «(z c , ( c , R CT ) .(46) 

, , J 8n^z c dt Mc 



a,b,c,d 

Here, it is understood that when the parton is a gluon, we need to replace f± and g\ with the momentum and helicity 
gluon distributions, respectively, as well as to have D\ — > G. The elementary cross sections are [42[ 

d! = 16ns 2 4 ^ M Xc ,x d ;x a ,x b ^xa,xb-,XaXd ' 

(allx's) 

rfAo- afc T_ >c T d _ 1 1 ,V ,">* / 47 x 

~ 167rs 2 4 ^ iw Xc^ d ;x a! Xi, iw xa,-xi,;-xc,x<i ' v*'.* 

(allx's) 

In the polarized case, a parton b is transversely polarized in a direction forming an azimuthal angle (f>s b around Pb 
and the transverse polarization of parton c forms the same azimuthal angle <ps r — 4>s b around Pc {4>s b and </>s c are 
defined analogously to (f>s B and 4> Rc , respectively, see Fig. 

Using Eq. (|45)l . the extraction of hi is possible either by using models for DiFF [Z(|, or by awaiting for e + e~ 
measurements to determine DiFF at BELLE [70j . in combination with the recent one in SIDIS at HERMES [37l ]. 
Experimental results at different energies can be presently related through DGLAP equations at NLL accuracy 
including full dependence upon R T [4l|. However, DiFF can be measured independently in the very same proton- 
proton collisions analyzed so far, by simply detecting another hadron pair in the other recoiling back-to-back jet. In 
fact, by generalizing the formalism of Eqs. (|4"Tj) and (l4"rH) , we have f42j 

dauu = A + cos (4> Rc - <j> Rn ) B + cos {2<j) Rc - 2cj> Rn ) C , (48) 
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where 



and 



6 = E / g| A»(».) ^ A'fo) dT ^ Cc i^r) Iggl igfr , a , 

C = E / gf /?(«•) A'fr) rfA ^f TdT ^ <^ c , Co, fl&r) ^ 5G<(z d , C, B? OT ) , (49) 



,b,c,d 



dAa, 



^! y 4 r , » M* . v v , (50) 



167TS 2 4 

(all X 's) 

with Xb = x a e~ vc e~ riD by momentum conservation. 

All functions A,B and C are interesting. The first two contain pairs of the DiFF D\,H^, one for each 
hadron pair: measuring the symmetric and the cos {4>r c — (pRo) asymmetric parts of the cross section for the 
pp — ► {hih2)j e tc{h\h2)jet D X process, it allows the extraction of D\ and and, in turn, of the transversity hi 
from the asymmetry An of Eqs. (|45I46[) in the corresponding polarized collision pp^ — > [hih^jX, Finally, the C 
observable describes the fragmentation of two transversely (linearly) polarized gluons into two transversely (linearly) 
polarized spin-1 resonances, which would not be available in a SIDIS process using spin-i proton targets. 



VI. CONCLUSIONS 



Using energetic hadronic (pion) beams on a transversely polarized NH3 target, it opens a new window at COM- 
PASS for the exploration of the partonic (spin) structure of the nucleon. When releasing the condition for collinear 
elementary annihilations, the leading- twist cross section for the production of Drell-Yan muon pairs contains several 
azimuthal asymmetric terms, whose combined study would shed light on yet unresolved puzzles. For example, the 
cos 20 term in the unpolarized cross section involves twice the exotic partonic density h^- , the socalled Boer-Mulders 
function, that describes the unbalance in the distribution of tranversely polarized partons inside unpolarized spin-i 
hadrons. It could give a natural explanation of the role of this contribution in the yet unexplained violation of the 
socalled Lam- Tung sum rule. The same hi happens in the sin(0 + (j> s ) asymmetric part of the polarized cross section, 
with <p s the target polarization orientation. It is convoluted with the tranversity hi, then offering an alternative strat- 
egy for its extraction with respect to the standard Collins effect in SIDIS. Finally, the polarized cross section contains 
also a sm((j) — </> s ) asymmetric part where the usual momentum density /1 is convoluted with fyj., the Sivers function, 
that describes the distortion of unpolarized parton distributions by the transverse polarization of the parent hadron. 
Studying all these partonic densities allows to build a 3-dim. map of partons inside hadrons via the determination of 
their angular orbital momentum, contributing to the solution of the yet pending puzzle about the spin sum rule of 
the proton. Moreover, a theorem based on very general hypotheses, would predict a sign change for f^~ T and h± being 
extracted from Drell-Yan measurements with respect to a SIDIS one. An experimental verification of this represents 
a formidable test of QCD in its nonperturbative domain. 

Monte Carlo simulations of the above azimuthal (spin) asymmetries at COMPASS, reveal that a very good statistics 
can be reached in short running time, given the relative abundance of energetic pions in the beam and the foreseen high 
luminosity. As a consequence, COMPASS is an ideal place for studying the sin(0— (j) s ) single-spin asymmetry involving 
the Sivers function, and for testing its predicted sign change with respect to SIDIS. Very small statistical error bars 
are reachable also for the other two asymmetries, but experimental constraints coming from older experiments (NA10, 
E615,..) reduce the size of the simulated asymmetry; a lack of parametrizations for the Boer-Mulders function, cither 
for proton targets or for pion beams, further makes its extraction an issue under debate. 

We have also shown that in semi-inclusive (transversely polarized) collisions with production of one hadron pair in 
the same jet, it is possible to isolate the convolution fi®hi® H^, involving the describing the fragmentation of 
a transversely polarized parton in the two observed hadrons, through the measurement of the asymmetry of the cross 
section in the azimuthal orientation of the pair around its cm. momentum. In the production of two hadron pairs 
in two separate jets in unpolarized collisions, it is possible to isolate the convolution fi ® fi ® ® H^, through 
the measurement of the asymmetry of the cross section in the relative azimuthal orientation of the two pairs. Since 
no distribution and fragmentation functions with an explicit transverse-momentum dependence are required, there 
is no need to consider suppressed contributions in the elementary cross sections included in the convolutions and 
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the discussed asymmetries remain at leading- twist. Therefore, contrary to what happens in SIDIS, proton-proton 
collisions offer a unique possibility to determine self-consistcntly all the unknown parton densities that are necessary 
to extract the transversity h\. We believe that this option should be seriously considered at COMPASS. 
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